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The double beta decays of 100 Mo and 150 Nd were studied 
in a Time Projection Chamber located 72 m underground. A 
3275 h exposure of a 16.7 g sample of metallic Mo enriched 
to 97.4% in 100 Mo resulted in a two-neutrino half-life, Ty 2 = 
(6.82±o53±0.68)xl0 18 y. Similarly, a 6287 h exposure of 15.5 g 
of Nd 2 3 enriched to 91% in 150 Nd yielded Ty 2 = (6.75t£g ± 
0.68) xlO 18 y. Lower limits on half-lives for neutrinoless decay 
with and without majoron emission also have been measured. 

PACS numbers: 21.10.Tg, 23.40.-s, 27.60. +j, 27.70. +q 



I. INTRODUCTION 



Double beta decay, the single isobaric jump of two 
units in atomic number, 



(A,Z) — ► (A, Z + 2) + 2e" 



■2z/ n 



is a rare second-order weak transition, directly mea- 
surable under favorable circumstances by means of the 
ejected pair of (3 particles. The phenomenon has been 
observed, with its two accompanying neutrinos (/3j32u), 
in several nuclei for which single beta decay is strongly 
inhibited or energetically forbidden |l]]. The search for a 
variation lacking the usual pair of neutrinos in the final 
state, is underway in several laboratories as a uniquely 
sensitive probe of the mass and charge conjugation prop- 
erties of the neutrino. Should a neutrinoless double beta 
decay {(3(3q v ) branch be seen, it would demonstrate that 
at least one neutrino is a massive Majorana particle 
at variance with the lepton conserving, massless Dirac 
particles of the standard model. 

A second nonstandard mode sometimes considered is 
neutrinoless decay with majoron emission ((3(3ov iX ). The 
majoron is a hypothetical boson, coupling to the neu- 
trino with sufficient strength in recent models, to make 
a significant contribution to the /3/3-decay rate |j-rcj. 

The three decay modes are distinguishable experimen- 
tally by the spectrum of the sum of the two /3-particle 
energies. The (3(3q v mode is characterized by a distinc- 
tive line spectrum at Q/3/3, whereas the flfov electron sum 
spectrum is a broad distribution, peaking at about 1/3 of 
the value. The various proposed P/3q v<x spectra are 
also broad distributions, generally distinct from fifov 

While all (3(3 decays observed by direct counting ex- 
periments have had the features of standard /3/?2i/> strin- 
gent lower limits on (3(3§ v and (3(3ov, x half-lives have been 



achieved. Geochemical and radiochemical (3(3 experi- 
ments, and searches for 7-rays following j3/3 decay to ex- 
cited levels of the daughter, do not distinguish among 
decay modes, but when their rates have been compared 
with direct-counting limits on the exotic modes, a pre- 
dominance of (3/32 V has always been implied [Q. 

The extraction of upper limits on neutrino mass and 
neutrino-majoron coupling strength (or of actual values 
for these parameters, should the /3/?ot/ and (3(3ov, x half- 
lives ever be measured) depends on complicated nuclear 
matrix element calculations. Comparison of theoretical 
and experimental half-lives for the observed fifiiv mode 
is an important test for these calculations Q]. 

Here we describe our final measurements |l of two 
relatively fast (3(3 emitters, 100 Mo (Qgg = 3034 ± 6 keV) 
and 150 Nd (Qpfs = 3367.1 ± 2.2 keV) g. We also discuss 
the disappearance of a persistent excess of high-energy 
/3/3oy,x"hke events that were present in our earlier data. 



II. EXPERIMENTAL APPARATUS 
A. The TPC 

The measurements were carried out in a time projec- 
tion chamber (TPC) located in an underground valve 
house in a canyon wall at the Hoover Dam. This site pro- 
vided a minimum of 72 m of rock shielding, and reduced 
the cosmic ray muon flux by a factor of approximately 
130 from that at the Earth's surface. 

The TPC has been described in detail in Refs. jb||ll). 
It is a rectangular polycarbonate (Lexan) box of inside 
dimensions 88 cm x 88 cm x 23 cm. The (3(3 source plane 
bisects the volume into two 10-cm-deep drift regions, and 
serves as the central drift field electrode. Planes of sens- 
ing wires parallel to the source are located near the TPC 
walls. Anode and cathode wires are mutually perpen- 
dicular, and provide x and y coordinates. The z coordi- 
nate is provided by the arrival time of drifting ionization 
electrons at the anode. Spatial resolution in x and y is 
determined by the wire spacing, 5.1 mm. The 5.0 mm 
z-resolution is established by the 5 mm/^s drift velocity 
and 1 MHz frequency of wire read-out. The correspond- 
ing 1 fis "time buckets" are clocked into an 80-deep shift 
register, in which the full 10 cm drift distance spans the 
first 20 buckets. The remaining 60 buckets are used to 
record any hits which occur in the 1 ms following the ini- 
tial trigger. This scheme permits tagging background 
from 214 Bi events by recording the subsequent 164-/is 



1 



pp 



20 cm 



Side A 



Side 



B 



TPC Wall 



Source 



Grid 

Field Wires 
Anode Wires 
Cathode Wires 



FIG. 1. Schematic representation of the source plane and 
wire array geometry. The sinusoids represent the x-z view of 
helical tracks from two electrons emitted by the source. 



214 Po a-particle. The required drift field is established 
by a negative potential applied to the source plane. Drift 
field and anode wire potentials are compensated contin- 
uously for fluctuations in barometric pressure. The ge- 
ometry of the source and wires is illustrated in Fig. [l| 

The TPC is shielded on all sides by 15 cm of lead. 
Notoriously radioactive materials, such as circuit board, 
solder, ribbon cable, circuit components, and connectors, 
are all completely outside of the lead shield. Interleaved 
with the lead shielding are six cosmic-ray veto panels in a 
47T arrangement which identifies muon associated events 
at a rate of approximately 0.25 Hz. 

A magnetic field is applied perpendicular to the source 
plane by a pair of coils mounted in an iron flux return. 
An electron emitted from the source follows a helical tra- 
jectory about the z-axis. The momentum of the particle 
and its angle to the magnetic field are determined from 
the parameters of the helix fitted to its track. The 1200 
gauss field (approximately twice the strength in our ear- 
lier experiments without the flux return) is uniform to 
±0.8% over the active volume of the TPC. 

Helium mixed with 7% propane flows through the 
chamber at a rate of 57 1/hr. Mixture proportions and 
flow rate are both controlled by a mass-flow controller 
system. Prior to mixing, the helium is passed through a 
liquid-nitrogen cold trap, which freezes out all measur- 
able Rn produced by Ra within the gas supply cylinders. 

The trigger used in the present experiment requires a 
start hit and at least one hit in buckets 3-5. In addition to 
complete tracks, the trigger will capture tracks having a 
z component as short as three time buckets. This feature 
increases the efficiency at which short a-particle tracks 
are saved, while maintaining some protection from noisy 
wires. Here we refer to a-particles other than the 214 Po 
type already included in 1 ms sensitive period following a 
214 Bi trigger. An a-particle appearing within an hour or 
two at the site of a (3(3 candidate event is a good indicator 
that the event was background from one of the primordial 
decay chains. The trigger fires at about 0.6 Hz. 

A small dead time in the chamber arises primarily from 
the cosmic ray veto. Triggers are vetoed for a period 
of 30 fis following a hit in any of the veto panels. The 



resulting dead time is about 0.6%. The overall dead time 
of the system is (0.87 ± 0.05)%. 

The data acquisition system is able to direct the latch 
modules to disable the occasional noisy wire. Typically, 
the TPC is operated with about 1% of the wires dead or 
disabled. 

Routines for off-line selecti on o f one- and two-electron 
events are discussed in Sec. VII. Only those 2e~ pairs 



with electrons emerging from opposite sides of the source 
are analyzed. 



B. Isotope Preparation 

The molybdenum, enriched to 97.4% in 100 Mo, was 
purchased as M0O3 from a German supplier, and is of 
Russian origin. The oxide is not the optimum form for 
detection of a-particles associated with background pro- 
cesses. As an example, the 7.7 MeV 214 Po a-particle 
range is diminished from about 22 mg/cm 2 in the metal 
to about 10 mg/cm 2 of Mo in the oxide. For this reason, 
the M0O3 was reduced to the metal in our laboratory in a 
quartz tube in a high temperature hydrogen atmosphere. 
About 20.0 g of fine metallic powder were available for 
source fabrication. 

The 150 Nd consisted of 3 g of Nd 2 03 purchased from 
the same German supplier, plus ra 1/4 of an 82 g batch 
of Russian Nd203 obtained through an agreement with 
the late Professor A. A. Pomansky of the Institute for 
Nuclear Research in Moscow. Elemental Nd is reactive, 
so it was left as the oxide. 

Samples of both isotopes were forwarded to Pacific 
Northwest Laboratories (PNL) for photon counting and 
mass spectroscopy (ICPMS). The PNL results indicated 
radioisotopes of europium in the 150 Nd at levels too great 
for /3/3-decay studies. The INR and commercial batches 
of neodymium oxide then were sent to the Ames Lab- 
oratory Materials Preparation Center, where they were 
combined and processed with a liquid chromatography 
system. The 150 Nd isotopic enrichment of the combined 
batches was 90.9%. Of the initial 85 g of Nd 2 3 , 94% 
survived the purification process. Samples were again 
forwarded to PNL where post-purification measurements 
indicated a reduction of more than two orders of magni- 
tude in europium contamination ]l2] | . This reduction was 
more than sufficient, although we were requested not to 
reveal the level of sensitivity of the measurements. Ura- 
nium and Thorium were not detected. 

Meanwhile our portion of the Nd203, which had been 
stored in sealed plastic bags, changed color from its char- 
acteristic robin's-egg blue, to a blue-grey. This indicated 
that some of the oxide had been hydrated to Nd(OH)3. 
This compound is undesirable for the same reasons dis- 
cussed in the case of the 100 Mo oxide. The hydrate was 
transferred to a quartz crucible and heated to 1100 °C in 
an electric furnace. This restored the characteristic blue 
color, and the approximately 12% loss in mass observed 
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is consistent with conversion of Nd(OH)3 to NC12O3. Af- 
ter processing, approximately 19.0 g remained. It was 
stored in an evacuated jar until needed. 



C. Source Preparation 

A thin layer of 100 Mo was deposited on each of two 
4 /zm aluminized polyester (Mylar) substrates, each sup- 
ported by a half-section of a Lexan source frame assem- 
bly. A nitrogen gas powered gun was used to inject 
the powdered isotope material at high speed into a glass 
box positioned in turn over each substrate. The result- 
ing cloud of powder settled onto the Mylar, leaving a 
thin, uniform deposit of isotope. This fragile layer of 
powder was then fixed to the Mylar with a fine misting 
of polyvinyl acetal resin (Formvar) solution. A 1 pCi 
droplet of 207 Bi solution was applied to the geometrical 
center of one of the source halves in a 3 cm square region 
left devoid of powder by a mask. After drying, the Bi was 
fixed in place with a drop of Formvar solution to prevent 
migration to the (3(3 deposit. (The conversion electrons 
emitted by the 207 Bi provide a continuous monitor of de- 
tector performance.) When dry, the two substrates were 
placed face-to-face to form a sandwich of Mylar and iso- 
tope, with the isotope in the interior, and the aluminized 
surface of the Mylar facing outward. 

All seams between the Mylar and the Lexan frame were 
sealed with low activity epoxy, isolating the interior of the 
source from the TPC gas. The interior of the source was 
vented to atmosphere through an oil bubbler, allowing 
the slight overpressure of the TPC to press the sheets of 
source Mylar tightly together into a thin, flat plane. Both 
the 100 Mo and 150 Nd (3(3 sources were prepared in the 
same manner. The masses of the isotope deposits on the 
two sources were (16.7±0.1) g of Mo and (15.5 ±0.1) g of 
Nd2C>3. Full details of the source preparation procedure 
may be found in Ref. @]. 



III. DETECTOR PERFORMANCE 

A. a-Particle Detection Probability 

Detection of a-particles is an important aid in identifi- 
cation of /3/3-decay backgrounds. The thicknesses of the 
two isotope deposits were chosen to assure that a large 
fraction of a-particles produced on or in the (3(3 source 
would escape to the chamber gas where they would be ob- 
served. The probability P Q that an alpha particle from 
the source will enter the TPC gas has been calculated 
with the CERN Library GEANT Monte Carlo, applied 
to the geometry and composition of the source. P Q was 
also measured, for the case of the 214 Po a particle, fol- 
lowing an injection of 222 Rn of sufficient strength to over- 
whelm pre-existing activity. Here P Q ( 214 Po) was taken 
as the number of single electrons of energy > 1.5 MeV 



(essentially all were 214 Bi (3 particles) that had an alpha 
track in the following millisecond, divided by the sum 
with or without an alpha. For details, see Ref. ||. The 
214 Bi from the Rn injection settles on the source surface, 
but as long as the source thickness is less than the alpha 
particle range, P Q determined in this way is a good ap- 
proximation for 214 Bi located anywhere on or within the 
source. The results of these methods are given in Table @. 

To the extent that 214 Bi is on the surface of the source 
plane (»100 % for the Mo source and «50% for the Nd 
source, Sec. VIA), the 214 Po a-decay can be identified 



with higher efficiency than implied by these escape proba- 
bilities. For surface 214 Bi, even the a-particle that buries 
itself completely in the source, is nearly always accom- 
panied by the release of "shake-off" electrons, which pro- 
duce a one- or two-bucket "blip" at the x,y location of the 
decay. Observation of such a blip in the 1 ms post-trigger 
interval identifies the trigger event as a 214 Bi decay. This 
association is confirmed by the delay-time distribution of 
blips, which shows the characteristic 164 /is half-life of 
214 Po. 

Since the blips substantially increase the identification 
chances for the 214 Po a particle (see Table |) they are 
used as well as a-particle tracks in the routine rejection of 
214 Bi. The alpha detection enhancement from shake-off 
electrons on the 150 Nd source is smaller because of sub- 
surface 226 Ra present in the 150 Nd2O3. Other a-decays 
used for background estimation are 212 Bi and 212 Po, also 
included in Table 1. 



B. Opening Angle Acceptance 

The TPC's inability to reconstruct tracks making an- 
gles to the magnetic field with cosines near or 1, the 
requirement that electrons from a /3/3-decay emerge from 
opposite sides of the source, and multiple scattering of 
the electrons within the source, all conspire to distort 
the measured opening angle distribution and obliterate 
the distinction between opening angle distributions for 
(3(3 decay and background. Because attempts to salvage 
the distributions through Monte Carlo corrections have 
not been very successful, we do not make use of opening 
angle information in our present analysis. 

However, it is worth pointing out that the requirement 
for opposite-side electrons leads to low efficiency for small 
opening angles, and thereby helps suppress one of the 
potential /3/3-decay backgrounds, Moller scattering (see 
Sec. VI). Moller opening angles are in the range cos$m ~ 
+0.2 to +0.5. The efficiency for reconstructing events in 
this range is relatively low, and will selectively suppress 
Moller events. 
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TABLE I. Ranges and escape probabilities for several important a-particles. Values in parenthesis for 214 Po refer to a-decay 
identification probability when "blips" from shake-off electrons are used in addition to full a-particle tracks. Energies are taken 
from Kaplan |13||. Measured P a values from Ref. |J have been corrected for loss of 214 Bi events when the a comes within ~ 15 
(is and spoils the /3 track, or comes after the 1024 fj,s window. 



a-Decay 


a Energy 
(MeV) 


Mean Range 
in air (cm) 


100 

Measured 


Estimated Escape 
Mo source 

GEANT 


Probability, P a (%) 

150 Nd source 
Measured GEANT 


214 Po 


7.6804 ± 0.0009 


6.907 ± 0.006 


73 ±5 


75 


72 ±5 


76 


212 Bi 






(98 ± 2) 




(85 ± 7) 




6.0466 ± 0.0027 


4.730 ± 0.008 


60 ±7 a 


65 


58 ±7 a 


66 


212p 


8.7801 ± 0.0040 


8.570 ±0.007 


79 ±4 a 


79 


78 ±4 a 


80 



a Scaled from the Po measurement by the Bragg-Kleeman approximation for a-particle ranges in materials. 
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FIG. 2. The measured energy spectrum of 7x 10 5 207 Bi IC 
electrons. Inset figure shows the weak peak near 1.7 MeV. 
Electrons with tracks nearly parallel to the magnetic field 
lines ( | cos i9 1 > 0.8) have been cut. Conversion line energies 
and relative strengths are represented by vertical lines. 



C. Energy Resolution 

A convenient line source for displaying the energy res- 
olution of the TPC is 207 Bi, which provides internal 
conversion (IC) energies at approximately 0.5, 1.0, and 
1.7 MeV. Although the dozen conversion electrons per 
hour emitted from the 1 pCi 207 Bi deposit in the center 
of the /3/3-decay source is sufficient for day-to-day mon- 
itoring of detector performance, a stronger source was 
required temporarily for accurate energy resolution stud- 
ies. This is particularly true for high energies, since the 
1.7 MeV 207 Bi conversion line is extremely weak. There- 
fore, a separate source assembly was built and peppered 
with 16 deposits of 207 Bi totalling 6 nCi, distributed over 
the area normally covered with powdered /3/?-decay iso- 
tope. 
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FIG. 3. Energy resolution, in percent, of the TPC as a 
function of electron energy and the cosine of the angle that 
electron makes with the magnetic field lines. This function is 
derived from studying the conversion lines of 207 Bi. 



The resulting energy spectrum of 7xl0 5 events is shown 
in Fig. ||. Energies of electrons with tracks nearly par- 
allel to the magnetic field (|cosi?| > 0.8) are less well 
determined, and have been excluded from this spectrum. 
The widths of the peaks are determined by fitting each of 
them with a pair of Gaussians, centered on the K and L 
lines. The corresponding a is 26 keV, 39 keV and 68 keV 
for the 482 keV, 976 keV and 1682 keV K lines, respec- 
tively. The width of each conversion line can be measured 
at various cosines, and a model of resolution versus co- 
sine and energy can be constructed. Figure ^ shows the 
results of this analysis. This function provides a statis- 
tical estimate of the uncertainty in the energy of a given 
electron track. The Monte Carlo program, described in 



Sec. HIE, uses this function to randomly smear the en- 
ergies of the simulated events. 
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D. Energy Threshold 

The energy of a particle is determined by fitting a helix 
to the three-dimensional reconstruction of its trajectory. 
Low-momentum particles will create helical tracks of rel- 
atively small radius and low pitch. As the diameter of 
the helix approaches the x- and y-resolution of the TPC, 
the quality of the reconstruction diminishes. Similarly, 
if the z-separation of adjacent cycles of the helix is less 
than the ^-resolution of the TPC, the helix will not be 
resolved, and only a single, wide swath of ionization will 
be reconstructed. 

Although some sensitivity clearly remains for energies 
as low as the 57 keV 207 Bi Auger electron jll]], for this 
analysis a 250 keV energy threshold was imposed. The 
TPC does not reach optimum efficiency until somewhat 
higher energies, but in place of still higher energy cuts, 
which would sacrifice a significant amount of good data, a 
model of the detector response was determined and used 
in subsequent analyses. 

Also included in the model is the high-energy response 
which is limited as well, because the pitch and radius of 
the track tend to increase with higher energies. At some 
point, the track becomes too straight for the fitter to 
determine reliable helix parameters. Fortunately, a large 
fraction of events satisfy the requirements for good fits, 
and can be very nicely reconstructed. 



E. Detector Response Model 

One could determine the detector response to a partic- 
ular (3 spectrum by placing a calibrated (3 source in the 
TPC, and noting the ratio of the measured spectrum to 
the true spectrum as a function of energy. Such a func- 
tion would not be generally applicable, however, because 
its dependence on the energy derivatives cannot be ex- 
tracted from a single test spectrum. For this reason, it 
is useful to separate the detector effects into two classes 
— intrinsic and instrumental, and deal with the former 
by Monte Carlo. The intrinsic effects include scatter- 
ing and energy loss of electrons within the (3(3 source, 
and backscattering of electrons from the TPC walls. In- 
strumental effects begin with the collection of ionization 
electrons from the gas, and include such factors as the 
sampling cell size, electron attachment by impurities, gas 
amplification, the readout electronics, and track recon- 
struction software. 

Intrinsic effects were modeled with the GEANT code. 
The model could then be applied generally for any (3 
spectrum. More importantly, the model can predict tea- 
electron behavior where the electrons are correlated in 
energy and opening angle, as for example (3(3 decay and 
Mollcr scattering. 

Instrumental effects are too complicated for a reliable 
Monte Carlo. Their two consequences are the resolution 
function, already determined from 207 Bi measurements in 
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FIG. 4. The normalized detector response as a function 
of energy and cos(i9), the angle between the track and the 
magnetic field. This is the electron survival probability after 
its emergence from the source plane. 



Sec. [II C , and the probability that an electron that has 
already survived the GEANT process for intrinsic effects 
will also survive the track reconstruction process. Since 
this probability is the same for any electron emerging 
from the source with a given energy and polar angle, it 
can be deduced from a comparison of a TPC-measured (3 
spectrum with the resolution-smeared output of the cor- 
responding intrinsic-effects Monte Carlo. This was done 
with the equilibrium 90 Sr, 90 Y combination (3 spectrum. 

Input electrons for the Monte Carlo were sampled from 
the 90 Sr and 90 Y spectra of Ref. @. The output en- 
ergies then were smeared by the resolution function of 
Sec. [IIC. Comparison with the spectrum of recon- 
structed electrons from a calibrated drop of 90 Sr solu- 
tion applied to a dummy (3(3 source in the TPC, revealed 
the survival probability of an electron after emerging 
from the source plane (Fig. ^|.) This survival probability, 
extrapolated for energies above 1.5 MeV, was then in- 
cluded in the GEANT Monte Carlo and used to regener- 
ate the 90 Sr lone electron spectrum; the results compared 
with the corresponding measured spectrum are shown in 
Fig. | 

Measured 2e events, generated by Moller scattering 
in the source plane, can also be compared with the sim- 
ulated results. The ratios of scanned 2e~ events to lone 
electron events, both taken from measured 90 Sr decays, 
were calculated at nine singles thresholds. As seen in 
Fig. ||, the measurements are in agreement with the 
Monte Carlo. Also, from this figure, a systematic un- 
certainty of 10 % can be estimated and is applied to all 
Monte Carlo results involving electrons. 
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FIG. 5. The measured 90 Sr lone electron spectrum (his- 
togram) and the simulated spectrum (solid line.) The sim- 
ulation includes energy smearing and the detector response 
model. 
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FIG. 6. The ratio of 2e~ Moller events to lone electron 
events as a function of singles threshold. Polar angles (•&) are 
restricted to | cosi9| < 0.8. The measured data are shifted by 
-0.01 MeV for clarity. 



IV. TESTS OF DETECTOR RESPONSE MODEL 

The response model of the previous section can be used 
to calculate several specific detection efficiencies which 
can be tested by direct measurement, as summarized in 
Table [n[ 

A. Single-Electron Efficiency, e e 

1. e e from 207 Bi Internal Conversion 



electrons from the two lines and the corresponding num- 
ber of reconstructed tracks. An independent measure- 
ment was made with the much hotter 6 nCi 207 Bi source 



described in Sec. IIIC . The activity of this uncalibrated 
207 Bi source was determined from the TPC trigger rate, 
which was dominated (98%) by the 207 Bi. Corrections 
were included for X-rays, Auger electrons, etc. This 
method yielded e e = (48.6 ± 0.8)% for the 1 MeV line. 



2. e e from 214 Bi /3-decay 

A TPC measurement of the 214 Bi (3 spectrum (Fig. 0) 
provides another efficiency check. Approximately 1 nCi 
of 222 Rn was injected into the TPC to produce a de- 
posit of daughter products, including 214 Bi, on the source 
plane. Since each 214 Bi /3 decay is followed by the a de- 
cay of 164 /is 214 Po, the requirement of a spatially cor- 
related a particle in the millisecond following the trigger 
was used to select triggering 214 Bi events, without re- 
gard to the quality of the (3 tracks. From the resulting 
46 15 1 eve nt sample, the single electron analysis software 
(Sec. VII ) reconstructed 9561 events originating on the 
source plane with energies above 500 keV. The energy 
cut was imposed to avoid complications related to the 
threshold of the detector and to a 5.5% contamination 
of the sample by 214 Pb events. A correction factor of 
1.030 must be applied to the 9561 events to account for 
(3 tracks spoiled by overlapping a tracks that satisfy the 
a selection criteria, but have short delays. 

Only 51.9% of the true 214 Bi (3 spectrum survives a 
500 keV cut. However, the primary sample was already 
depleted in low energy events, since a Monte Carlo cal- 
culation predicts that (8.1 ± 0.8)% of the electrons fail 
to escape the source and cause a trigger. An additional 
28.2 ± 2.3% correction eliminates contamination of the 
sample by (3 tracks originating on the TPC walls rather 
than the source plane. The resulting single-electron effi- 
ciency e e = 



9561 x 1.030 x [1 - (.081 ± .008)] x 100 
46151 x .519 x [1 - (.282 ± .023)] 



= (52.6±1.6)%. 



The calibrated 207 Bi source, at the center of the 150 Nd 
isotope deposit, has an activity of 140.0 ± 4.0 decays per 
hour. Given the conversion rate of the 1064 keV 7-ray, 
a 0.6% loss due to electron tracks spoiled by X-rays and 
Auger electrons, and a further 2.2% loss due to an ac- 
companying second conversion electron, 13.4 ± 0.4 sin- 
gle K,L,M... electrons from the 1064 keV transition can 
be expected every hour. The same calculation for the 
570 keV 7-ray gives 2.98 ± 0.09 single electrons per hour. 
During 6286.6 hours of livetime, 39 724 electrons were ob- 
served in the energy range expected from conversion of 
the 1064 keV 7-ray, and 8907 electrons for the 570 keV 
transition. Thus an overall efficiency of e e =(47.2± 1.2)% 
can be calculated by combining the number of expected 



B. Two-Electron Efficiency, e ee 

1- See from the 207 Bi decay scheme 

The small 207 Bi deposit located in the center of each 
source assembly can also be used to determine the 2e _ 
efficiency. Occasionally, two transitions in the cascade 
will each yield a conversion electron, producing a 2e _ 
event. The product of the conversion rates for the 0.5 and 
1 MeV lines gives the rate of IC-IC pairs as 2.914 x 10 -3 
per 1 MeV 7-ray. This implies that an IC-IC cascade 
should occur at a rate of 2.157x10" 

-3 per 207 Bi decay The 
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FIG. 7. The TPC-measured 214 Bi /3-decay energy spec- 
trum following 222 Rn injection. Solid curve is the GEANT 
Monte Carlo prediction normalized above 0.5 MeV. The dif- 
ference arises largely from contamination of the event sample 
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FIG. 9. Measured 214 Bi 2e" spectra, satisfying a 0.2 MeV 
singles threshold, following 222 Rn injection. Gaussian 
ideogram (solid curve) is a better representation of the ex- 
pected signal than is the histogram. Some 214 Pb contami- 
nation is present, but is not included in the GEANT Monte 
Carlo predictions (crosses). 
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FIG. 8. Scatter plot of electron energies for events in the 
150 Nd 2e~ data set which emerged near the 207 Bi deposit. 
The circle is centered on 481.7 keV and 975.6 keV (the K con- 
version energies). Moller and Compton scattering contribute 
events outside the circle. 



activity of the drop on the 150 Nd source is 1.05±0.03 pCi, 
which should give 0.3017 ± 0.009 IC-IC pairs per hour. 

The 2e~ stripper used to extract /3/3-decay events also 
saves 207 Bi IC-IC decays which are easily tagged by their 
location at the center the source plane. In the scatter 
plot of Fig. ||, most of the points are IC-IC events, as 
is evident from the clustering near the 481.7 keV and 
975.6 keV conversion energies. These events are selected 
by a circular cut which includes 99% of the IC-IC signal. 
Points outside the circle [[llj represent Moller scattering 
of single conversion electrons, or single conversion elec- 
trons accompanied by a Compton electron from a cascade 
gamma ray. 

The selection method described above identified 253 
IC-IC events in a livetime of 6176.6 hours, giving an ob- 
served rate of 0.041 ± 0.003 events per hour. When com- 
pared to the estimated activity of the 207 Bi deposit, this 
rate implies an efficiency of e ee = (13.6±0.9)%. 



2. See from the 214 Bi decay scheme 



The 222 Rn injection described in Sec. IV A 2 yielded 
a sample of 46 151 214 Bi decays from which the offline 
analysis reconstructed a total of 189 2e~ events (Fig. |^). 
A single-electron threshold of 400 keV then was imposed 
on the data to escape the influence of 214 Pb contamina- 
tion. The 64 above-threshold events were individually 
examined, and 59 good 2e~ events were identified. 

These 59 events included contributions from (3 + IC, j3 
+ Moller, (3 + Compton, etc. (See Table [n]). Unlike the 
situation in Fig. |^, however, the events involving IC are 
not easily isolated and counted. Therefore, a comparison 
of the j3 + IC subset to the number of (3 + IC events 
expected from the decay scheme is not a viable method 
of calculating the efficiency. Instead, one must take the 
quotient of all 59 measured 2e~ events and the total of all 
types of 2e~ events produced in the source. An estimate 
of the latter requires the very response model that we 
are trying to test by this measurement. To avoid this 
circular argument, the 214 Bi 2e~ efficiency is rejected as 
a useful test. 

Nevertheless, one can still carry out a meaningful check 
of the model by directly comparing the 59 measured 2e~ 
events with the number that the model predicts would be 
measured from 46151 2 Bi decays. The overall conver- 
sion rate was calculated by incorporating the /3, 7, and 
IC characteristics of the complicated 214 Bi decay scheme 
in a stand-alone Monte Carlo program. This simulation 
considered >99% of the allowed decay channels, and es- 
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TABLE II. TPC-measured efficiencies and two-electron 
events as tests of the Monte Carlo model. 





Single electron 


Two-electron 




efficiency e e (%) 


efficiency e ee (%) 




Measured Monte Carlo 


Measured Monte Carlo 


207 Bi 


53.0 ±0.1 


13.9 ±0.1 


1 pCi 


47.2 ± 1.2 


13.6 ±0.9 


6 nCi 


48.6 ±0.8 


Two-electron events 


214 Bj 


52.6 ± 1.6 54.6 ±0.2 


59 ± 8 52.9 ± 3.4 



timated the conversion rate to be 0.0194 single conver- 
sion electron per 214 Bi /3-decay (Table IV). Feeding this 
information into the response model, with the usual cor- 
rections for wall events and spoiled tracks described in 
Sec. 1VA2, resulted in a prediction of 52.9 ± 3.4 2e~ 
events from the 46 151 decays. The favorable compari- 



son to the 59 measured events is included in Table O. 



C. The Model Testing Process 

We test the reliability of the Monte Carlo for efficiency 
and 2e~ event prediction in those specific cases that have 
been measured directly, namely 207 Bi and 214 Bi as de- 
scribed above. The angular distributions and energies 
of fifov are replaced in the Monte Carlo with those for 
the appropriate isotope. For example, the opening angle 
distribution for the IC-IC events of 207 Bi [jl5| , 

dw/d(costf) - 1 ± 0.271 cos 2 d 

and the energies of the two conversion electrons, when in- 
serted in the Monte Carlo, give e ee for 207 Bi as 13.9 ±0.1, 
vs. a measured value of 13.6 ±0.9. The Monte Carlo and 
measured values for e e and e ee are summarized in Ta- 
ble ||. These results are consistent with the 10% sy stem- 
atic uncertainty assigned to the Monte Carlo in Sec. Ill E . 



V. THE (3fh v EFFICIENCY, epp 

Determination of the 100 Mo and 150 Nd half-lives from 
the observed rates requires knowledge of the (3(32v effi- 
ciency, epp. Since there is no direct way to measure epp, 
it must be estimated by Monte Carlo. The code em- 



ploys the detector response model (Sec. Ill E) and utilizes 
the opening angle and energy distributions expected for 
f3p2 V in the Primakoff- Rosen approximation p6| . The 
Monte-Carlo-generated efficiencies, as a function of sin- 
gles threshold, are shown in Fig. ITol. 



VI. BACKGROUNDS 

Candidates for /3/3-decay are characterized by two neg- 
ative electrons with a common point of origin on the (3(3 
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FIG. 10. The efficiency epp as a function of singles thresh- 
old. The statistical uncertainties are negligible compared with 
the 10 % systematic uncertainty (not shown). 



source. Although this signature rejects the vast majority 
of events from radioactive impurities, there are certain 
non-/3/3 processes that mimic the 2e~ topology. These 
include Moller scattering of single electrons born in the 
source, and (3-"f cascades in which a 7-ray Compton scat- 
ters, produces a photoelectron, or internally converts. 
(See Table [II.) Pair production (e + e~) is not a back- 
ground since the helical track of the positron would be 
rejected by a pitch of the wrong sign. 

Dominant among these mechanisms is (3 decay ac- 
companied by an internal conversion (/?±IC), a rela- 
tively common phenomenon in the 238 U and 232 Th de- 
cay chains. These families always are present at some 
level in the (3(3 isotope stock. In addition, daughters of 



gaseous 220 Rn or 222 Rn can migrate to the (3(3 source sur- 
face from Ra impurities elsewhere in the chamber or in 
the gas-handling system. 

The decay chains supply the major radioactivity feed- 
ing the other background processes as well. Fortunately, 
each of these chains contains one or more isotopes with 
unique and identifiable decay characteristics, so good es- 
timates of total uranium and thorium contamination can 
be made from the observed activity. 



A. The Uranium Series 

The uranium chain includes several isotopes with high 
enough values to be of concern. A summary of their 
properties is included in Table |y|. 

Activity below 230 Th is quantified by the characteris- 
tic 164 fis (3-ot sequence from 214 Bi and 214 Po ("BiPo" 
events). The TPC's 1 ms trigger duration includes 99% 
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TABLE III. Breakdown of 2e~ events, from 10 /3-decays, satisfying a 0.25 MeV singles threshold. These were generated 
by the GEANT model of intrinsic effects using the simulated decay schemes. The columns denote combinations of 2e~ events 
consisting of the beta (/?), Internal Conversion (IC), Moller (M), Compton (C) and Photoelectric (P) electrons. * = /3 or IC or 
M or C. Combinations with fewer than 10 events are not shown but are included in "All". In columns lacking a /3, the /3 did 
not escape the source. 



100 Mo 150 Nd 



Isotope 


/3 + IC 


P + M 


p + c 


IC + IC 


P + * 


All 


/3 + IC 


P + M 


P + C 


IC + IC 


P + * 


All 


234 Pa m 


548 


260 


1 








811 


546 


255 


1 





1 


804 


214p b 


2568 


2 








24 


2594 


2693 


3 








38 


2734 


214 Bi 


2163 


206 


167 


14 


21 


2577 


2107 


164 


151 


14 


36 


2480 


228 Ac 


972 


73 


50 


14 


21 


1135 


995 


72 


59 


18 


46 


1195 


212p b 


10 














10 


6 














6 


212 Bi 


234 


251 


9 


6 


1 


503 


241 


231 


17 


4 


3 


497 


208 rpj 


4501 


143 


204 


37 


28 


4924 


4527 


113 


197 


44 


78 


4966 


60 Co 








6 








7 








6 





3 


10 



TABLE IV. Properties of isotopes in the U, Th decay 
chains capable of producing high-energy /3+IC events. "To- 
tal" includes events with one or more conversion electrons, 
and "Single" indicates one conversion electron only . These 
are the raw conversion rates derived from Monte Carlo simu- 
lations of the decay schemes, without reference to the TPC. 





Q- Value 


IC per /3-decay (%) 


Isotope 


(MeV) 


Total 


Single 


234p a m 


2.29 


1.27 


1.18 


214p b 


1.02 


40.6 


35.1 


214 m 


3.27 


1.95 


1.94 


228 Ac 


2.14 


80.3 


67.5 


212p b 


0.573 


41.5 


40.2 


212 m 


2.25 


0.343 


0.339 


208 rpj 


4.99 


7.07 


6.91 



of the 214 Po decays. Alpha tracks that extend into the 
TPC gas are recorded and stored along with the j3 tracks 
responsible for the triggers. 

The 100 Mo source showed a low level of 214 Bi events 
distributed over the entire source plane, irrespective of 
the region where the 100 Mo was located. This uniformity 
suggests that very little of the activity in the 100 Mo runs 
arose from contamination of the isotope material. In- 
creased TPC gas flow rates were accompanied by fewer 
observed 214 Bi decays. The flow dependence implied that 
the primary source of 214 Bi in the 100 Mo data set was 
222 Rn in the chamber gas. Data collected during the 
150 Nd runs indicate that the INCI2O3 itself harbored ad- 
ditional 214 Bi activity roughly equal to that from 222 Rn 
contamination of the chamber gas. 

Quantifying the 214 Bi contamination, however, does 
not constrain 234 Pa m . A breach of secular equilib- 
rium could allow a 234 Pa m signal without any 214 Bi. 
Since the Mo refinement and enrichment processes could 
throw these isotopes out of equilibrium, the possibility of 
234p a m background must be considered. 



1. 214 Bi Background 

214 Bi with its 3.27 MeV Q-value, is potentially the 
single most dangerous background in /3/3-decay exper- 
iments with isotopes like 100 Mo (Q = 3.03 MeV) and 
150 Nd (Q = 3.37 MeV). The 214 Bi 2e" sum-energy spec- 
trum collected during the 222 Rn-injection runs (Fig. ||) 
is not unlike the expected [3 feu signal. The a particles 
and shake-off electrons from decay of the 214 Po daugh- 
ter a re exp loited to tag 214 Bi events with good efficiency 
(Sec. Ill A ) . The number of tagged events together with 
the known tagging efficiency P Q ( 214 Po) give an accurate 
measure of the number of untagged 214 Bi events remain- 
ing in the data. 



^untagged ^tagged (1 + 



l-^P Q ( 214 Po) 

P a ( 214 Po) 



where £ = 0.081 and r\ = 0.987 are corrections for cor- 
ruption of Bi events by Po a's in the first 15 /is, or for Po 
decays beyond the 1024 /is window. We found that we 
could also tag 214 Bi events by one or more of the three 
previous links in the uranium chain, with about 50% ef- 
ficiency in the more difficult 150 Nd case. The number of 
214 Bi events remaining after the primary 214 Po tag, how- 
ever, was already so small (Table |y|) that this more te- 
dious procedure was not needed. Another tool for reduc- 
ing the 214 Bi contribution is a cut on the strong 1.3 MeV 
conversion line in the singles spectrum, but that tech- 
nique also was not needed here. 



2. 214 Pb Background 

Although it has only a 1.0 MeV Q-value, 214 Pb has a 
very high IC probability which makes it a major source 
of low-energy background. The 214 Pb /3+IC singles- and 
sum-energy spectra were measured with the same 222 Rn 
inje ction d ata used to measure the 214 Bi spectrum in 
Sec. [V A 2. In this case, however, only manually scanned 







TABLE V. Summary of 2e~ /3/3-decay background events. 
228 Ac, 212 Pb and 212 Bi co ntamination are estimated from 
208 T1 activity (see Sec. [VI Bj). Also listed are Moller scattering 



contributions as estimated from lone-electron energy spectra. 
"Cuts" include a 250 keV singles threshold and 4-wire radius 
of exclusion about the 207 Bi spot at the center. 





Back- 




Number of Events 


/3/3-decay 


ground 




Tagged 


Remain 


Isotope 


Type 


Found 


Passed Cuts 


in Data 


100 Mo 


214p b 


24 


4 


0.16 ± 2.09 




214 Bi 


9 


4 


0.12 ± 0.11 




228 Ac 


- 




8.8 ± 3.5 a 




212 Pb 






0.08 ± 0.03 a 




212 B j 






0.6 ± 0.2 a 




208rpj 


14 


8 


5.8 ± 2.6 




Moller 






20.5 ± 3.6 ab 


150 Nd 


214p b 


76 


25 


33.3 ±10.2 




214 B j 


61 


44 


9.0 ± 4.7 




228 Ac 






11.9 ± 4.4 a 




212 Pb 






0.06 ± 0.02 a 




212 Bi 






0.8 ± 0.3 a 




208rpj 


25 


10 


7.9 ± 3.3 




Moller 






9.2 ± 3.1 ab 



a Monte Carlo systematic uncertainty included. 

b Excludes (3 + Moller events already included for the above 

isotopes. 
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events with no accompanying a-particle or shake-off elec- 
tron were considered. This sample contains < 2% con- 
tamination of 214 Bi events. 

The 214 Pb/ 214 Bi ratio of 2e" events from Table [fn] was 
used to determine the number of 21 Pb events in the (3/3 
data set from the number of tagged plus untagged 214 Bi 
events discussed in the previous section. (As a check, 
this ratio, measured from the 222 Rn injection data at 
a 0.200 MeV singles threshold, yielded 2.51 ± 0.30, in 
agreement with the corresponding Monte Carlo value of 
2.82 ±0.07 ±0.40.) 

It is possible to identify individually and remove many 
of the 214 Pb events in the /3/3-decay data set by tagging 
them with the characteristic "BiPo" signature that fol- 
lows with a mean delay of 28.4 minutes. 214 Bi decays 
occur at a rate of only a few per hour over the entire 
source area, so the chance of an accidental spatial and 
temporal coincidence between a (3(3 event and a 214 Bi 
decay is very small. If a 1 cm 2 area of coincidence and 
a two-hour search window are assumed, the chance is 
only about 1 in 500 that an event might be wrongly re- 
jected. The search has been carried out on all of the 
/3/3 candidates from both sources. For each candidate, 
the software searches the raw data file for any nearby 
event which shows a delayed hit pattern indicative of an 
a-particle, and occurs within the two hour (six half-life) 
time period following the candidate. After human con- 
firmation of a valid "BiPo" tag (Fig. |flj) the 214 Pb event 
is removed from the data. Results are shown in Table [v|. 



FIG. 11. An example of Pb /3±IC event identification. 
Top event is a low-energy /3/3-decay candidate. Bottom event 
is a "BiPo" sequence emerging from the same point on the 
source 82 minutes later, indicating that the candidate is prob- 
ably a 214 Pb /3±IC event. Tick marks near the vertex are at 
the same position in each figure. Open circles are the 214 Po 
a track, 747 fjs delayed. 



3. 234 Pa m Background 

Likely disequilibrium with the rest of the uranium 
chain forces us to evaluate the level of 234 Pa m contamina- 
tion by searching for its decays directly. This isotope has 
only one strong conversion line at 695 keV, so any 2e _ 
contribution should populate the singles-energy spectrum 
at that energy. There is no obvious excess observed at 
695 keV in either the 100 Mo or 150 Nd 2e~ singles-energy 
spectrum, so any 234 Pa m contribution is small. 

The 234 Pa m events that do occur must be preceded by 
the low-energy 234 Th /3-particle within several 1.2 min 
half- lives. Background identification techniques used to 
tag 208 T1 2e~ events will also find this decay sequence 
with high efficiency (although it would probably be la- 
beled as a 208 T1 event). Therefore, the background con- 
tribution from 234 Pa m is assumed to be negligible. 
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B. The Thorium Series 

The potential of thorium daughters for serious (3(3 
background is evident in Fig. 12 which shows measured 
singles- and sum-energy spectra for equilibrium 212 Pb 
+ 212 Bi -I- 208 T1, obtained by injecting 220 Rn gas into 
the center of the TPC. Also shown are corresponding 
spectra generated by a Monte Carlo. Noticeable in the 
sum-energy spectrum is a significant population in the 
3 MeV region where one looks for /3/?oi/-decay in 100 Mo 
and 150 Nd. 

The presence of thorium chain contaminants can be 
recognized from observation of the 0.3 fis, mass-212 
"BiPo" sequence. Within the temporal resolution of the 
TPC, both the (3- and a-decays will appear to occur si- 
multaneously, and will both be recorded in the trigger 
hits of the event. This signature can be imitated by «1% 
of the 214 Bi decays in which the 164 /is a-decay of 214 Po 
occurs in the first few microseconds. A small correction 
for the mass-214 effect can be applied once the uranium 
chain activity is estimated, and does not interfere with 
using 212 Bi as a thorium tracer. 

Unlike the uranium chain contaminants, thorium ac- 
tivity is not distributed over the (3(3 source by Rn in 
the gas supply. Gas is introduced slowly into the TPC 
from the top. As demonstrated by a separate injection of 
220 Rn into the TPC gas feed line, the charged daughter 
ions of this short-lived radon collect near the point of en- 
try, and subsequent activity is localized to the top of the 
chamber. Since this pattern is not observed in normal 
data runs, the thorium activity is within the chamber, 
probably in the (3(3 source itself. 

The serious background isotopes of the thorium se- 
ries are separated by a maximum half- life of 1.9 years 
( 228 Th). Considering the elapsed time between process- 
ing and measurement of our isotopes, it is safe to say that 
the thorium chain from 228 Ac to 208 T1 is in approximate 
equilibrium. 



1. 20S Tl Background 

With a Q-value of almost 5 MeV and 0.5% of decays 
generating two electrons passing the standard cuts (Ta- 
ble [II), 208 T1 is another particularly dangerous back- 
ground. Fortunately, the highest-energy decays are al- 
most exclusively from internal conversion of the 2.61 MeV 
7-ray, so events near 3 MeV can be removed by a cut 
on single electrons near the 2.53 MeV conversion energy, 
without significantly impacting the rest of the spectrum. 

The use of the mass-212 "BiPo" count to quantify the 
activity of the Th series is not as straightforward as in 
the mass-214 case. The effectively instantaneous 212 Po 
a-particle is often lost within the 212 Bi (3 track. It is 
sometimes possible, however, to identify a /3/3-decay can- 
didate as a 208 T1 event by a back search for the a-decay 
which must precede it with a 3.1 minute half- life. For all 




2 2.5 3 3.5 
Singles Energy (MeV) 

FIG. 12. Measured 212 Pb + 212 Bi + 208 T1 equilibrium 2e" 
spectrum, satisfying a 0.2 MeV singles threshold, following 
220 Rn injection. Solid curve is the Gaussian ideogram. The 
corresponding GEANT Monte Carlo prediction (crosses) in- 
cludes the branching and loss of 212 Bi events spoiled by the 

212 P0 Q. 



/3/3-decay candidate events, the preceeding 20 minutes 
(six half-lives) is searched by software for any particle 
originating near the candidate's vertex. Particle identi- 
fication is left to the judgement of the person analyzing 
the data. 

The efficiency of this search procedure is the probabil- 
ity of seeing the 212 Bi a-decay, P a ( 212 Bi) from Table Q. 
(This search would be enhanced by the conversion elec- 
tron from the 40 keV level in 208 T1 which occurs with 70% 
of the 212 Bi a-decays, but because of its low energy this 
electron rarely reaches the TPC gas.) The corresponding 
208 T1 count is included in Table 0. 



2. 21 Bi Background 

Although is 2.25 MeV, 2e~ production per (3 decay 
by 212 Bi is only 10% of that by 208 T1 (Table |llj]) . Branch- 
ing increases the 212 Bi portion to 18%, but finally it is 
reduced to 4% by corruption from the essentially simulta- 
neous 212 Po a-particle, which appears in the gas 80% of 
the time (Table Q). This small background contribution 
of 212 Bi is also included in Table [v|. 



3. 212 Pb Background 

The number of 212 Pb 2e~ events can be determined 
from 208 T1 decays in a manner similar to that for 212 Bi 
above. The branching fractions and 2e~ probabilities 
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predict that 212 Pb events should be detected < 0.6% as 
frequently as similar 208 Tl decays. These are also in- 
cluded in Table 0. 

J h 22S Ac Background 

The complicated decay scheme of 228 Ac almost always 
produces an internal conversion, and will frequently pro- 
duce multiple conversions. Decays with three or more 
electrons appearing will not be accepted as /3/3-decay 
events, so this characteristic helps reduce the probability 
of 228 Ac contributing to the measured /3/3-decay spectra. 

Although the singles-energy Monte Carlo spectrum for 
228 Ac /3+IC events has a distinctive peak near 0.85 MeV, 
there is no hint of such a peak in the measured 2e _ spec- 
tra. 

We can estimate the 228 Ac contribution from the 208 Tl 
activity as described in the previous sections. Assuming 
that the thorium chain is in equilibrium, branching frac- 
tions and 2e~ probabilities predict that 228 Ac 2e~ events 
should be detected with ~ 65% the frequency of similar 
208 Tl decays. These are also included in Table [v|. 

C. Other Sources of Conversion Backgrounds 

Moe and Lowenthal (l7) have investigated the possibil- 
ity of /3+IC backgrounds produced by isotopes outside 
of the natural decay series. They considered nuclei pro- 
duced from cosmic ray activity, contaminants from atmo- 
spheric nuclear weapons test fallout, common man-made 
materials such as 60 Co, and long-lived non-series isotopes 
such as 40 K. In order to be considered a serious source of 
/3+IC background, an isotope must have significant inter- 
nal conversion probability and Q- value, a half-life greater 
than the order of six months or, alternately, a plausible 
continuous production mechanism. (We have never seen 
a /3/3 candidate rate decline over the duration of the ex- 
periment.) None of the nuclides investigated satisfied all 
of these requirements. 

VII. ANALYSIS OF THE DATA 

A. Off-Line Event Selection 

Each day the data acquisition system recorded about 
50 000 events, from which a (3(3 signal of 2 or 3 events 
had to be extracted. Automated event stripping routines 
applied a set of empirical criteria engineered to select 
events which could be /3/3-decay candidates. A similar set 
of routines was designed for identifying lone electrons. 



B. Event Scanning and Fitting 

1. Lone Electron Events 

Since 100 Mo and 150 Nd are stable against single (3- 
decay, all lone electrons spawned in the source are due to 
unwanted radioactivity. Producers of 2e~ backgrounds 
generally emit far more copious numbers of single elec- 
trons. Studying the lone electrons^ emitted from the 
source plane helps establish limits on the number of asso- 
ciated 2e~ background events. This can be accomplished 
by fitting the lone electron energy spectrum with compo- 
nent background spectra to estimate each component's 
activity. 

Figure [ll| shows the energy spectra of lone electrons, 
whose track vertices are between 10 and 53 wires from 
the 207 Bi source, collected from all but 518 hours of the 
100 Mo runs and from the entire set of 150 Nd runs. These 
spectra were fitted with various background lone elec- 
tron spectra generated by the GEANT Monte Carlo us- 
ing Ref. [ pT[ for input spectra. Background component 
spectra considered were 40 K, 60 Co, 137 Cs and daughters 
of the 238 U, 232 Th and 235 U decay chains. In addition, 
108 Ag was considered for the 100 Mo data while 152 Eu and 
154 Eu were included in the fit to the 150 Nd data. 

The fit to the 100 Mo spectrum was of limited success. 
Difficulty reproducing the shoulder near 1.3 MeV is sug- 
gestive of a contribution from an unidentified contami- 
nant. We do not consider this fit to be valid. Since the 
likely heavy elements have been included, the unknown 
contaminant is not expected to be a strong IC emitter. 



1 It was his last wish before his death in 1990, that Prof. 
D. Skobeltzin of INR Moscow, hoped someone would investi- 
gate the anomalous scattering of (3 particles observed in his 
early cloud chamber experiments. Suggesting that a He-filled 
TPC was the ideal detector, Prof. G. Zatsepin relayed this 
wish to us, and requested that we publish this footnote. Sko- 
beltzin's photographs (eg. Nature 137, 234 (1936)) showed 
(3 particles from 214 Bi suffering apparent inelastic collisions 
in which they lost as much as 90% of their energy. As Sko- 
beltzin described it in a book honoring the 60th Anniversary 
of the birth of Prof. S. Vavilov published by the Academy of 
Sciences of the USSR in 1952, the phenomenon occurred in 
5-10% of the (3 tracks within 20 cm of the source, and was 
consistent kinematically with a ~ 10~ 10 sec in-flight decay of 
an unstable particle. We have never observed this effect in 
the many beta tracks examined in the course of our (3(3 work, 
and would certainly have seen it at the 5-10% level. Never- 
theless, to be quantitative, we specifically examined the first 
20 cm of 230 high-energy (3 tracks from the decay of 214 Bi. 
We saw no event with a sudden, large loss of energy. This re- 
sult rules out more than a 1% effect in our gas mixture, with 
90% confidence. Without seeing Prof. Skobeltzin's original 
photographs, we can only speculate that the anomaly resulted 
from nonuniform illumination of his chamber. 
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FIG. 13. Histograms of the lone electron energy spectra 
from 100 Mo and 150 Nd runs, and their fitted (at >200 keV) 
background components. These are 238 U (excluding 234 Pa m ) 
(dashed curve), 232 Th (dotted curve), the sum of 40 K, 60 Co 
and 137 Cs (dashed-dotted curve), and the total spectrum 
(solid curve). 



Its dominant contribution to 2e— background should be 
Moller scattering, which is accounted for by the lon e elec- 
tron spectrum, independent of any fit (Sec. VII C| ). 

The fit of the 150 Nd spectrum, however, yielded a good 
match _to the measured data. The results are shown in 
Although 212 Bi is in equilibrium with the 
it was fitted as an independent pa- 



Table [VI. 
232 Th decay chain 



rameter in order to calculate the 212 Po cv-particle escape 
probability. This was determined to be 67 .5^1%, within 
1.64a of the values listed in Table |. 

In addition to the isotopes listed in this table, oth- 
ers were considered and rejected. Daughters of 235 U, 
specifically 211 Pb and 207 T1, were acceptable to the fit 
but at the expense of 40 K. Ultimately, 235 U background 
was rejected by the fits to the 2e~ spectra discussed in 
Sec. VTID2. Similarly 152 Eu was also acceptable in the 
lone electron fits but not in the 2e~ counterpart. On the 
other hand, 154 Eu was rejected by the lone electron fit. 



2. Two Electron Events 

The offline stripper selects negative electron pairs 
emerging from opposite sides of the source plane. We 
did not reconstruct same-side events because of added 
complexity in event selection and fitting, as well as in 
understanding the associated efficiencies. 

The stripper attempts to reject events which are not 
/3/3-decay candidates. When it cannot reject an event 
with certainty, it will save it. This system ensures that 



TABLE VI. Activities corresponding to the fit to the 150 Nd 
lone electron energy spectrum of Fig. [131. Also shown are total 
2e~ events surviving standard cuts, as calculated from the 
activity, and from Table [v] which is considered to be more 
reliable. 





Activity 


Total 2 


e events 




Isotope 


(uBa/e) 


From Fit a 


From Table 


234 Pa m 


< 18.7 b 


< 8.7 b 






214p b 


71 ^+14 c 

/ l.J_4.1 


68 5 +1-8 

OO.O _8.o 


58.3 ± 10.2 


214 Bi 


Q 


58 3 +SM 


53.0 ± 


4.7 


210 B j 


4 3 + ,yn 8 

1-J—O.O 


o 






228 Ac 


98.3±i|;? d 


4i.i in 


11.9 ± 


4.4 a 


212p b 


(1 


n oi +0.09 


0.06 ± 


0.02 a 


208 rpj 


(1 


61.6 tlli 


17.9 ± 


3.3 


212 ffi 


20.4+^ 


0.7 it? 


0.8 ± 


0.3 a 


40 K 


43.71SI 









60 Co 


613 t| 4 


O i +0.7 
^•t _0.7 






137 Cs 


239 +l° 










a Monte Carlo systematic uncertainty included. 

b Limit at 90 % confidence level. 

c 214 Pb and 214 Bi were fitted together. 

d 228 212p b and 208 T1 wgre fitted together 



essentially all recognizable /3/3-decay events are accepted, 
and efficiency of the offline stripper is not an issue. Stud- 
ies conducted during the development of the stripper 
showed that it did not reject any events which would 
not have been rejected by a human scanner. 

The offline stripper saved 2e~ "/3/3-decay candidate" 
events at rate of about 25-30 per day of running. These 
events were individually studied by a physicist. Most 
events were rejected immediately as a multi-tracked 
backscattered lone electron, an indecipherable swath of 
ionization, or some other two-sided event which was 
passed by the stripper's policy of accepting events that 
it cannot categorize. 

Events that have the correct /3/3-decay topology are fit- 
ted with two independent helices. The scanning software 
determines the energy and opening angle of the event 
based on the parameters of the fit in conjunction with 
the magnet current, drift field voltage and atmospheric 
pressure information recorded with the event. If the fit 
is acceptable, the person scanning the data can elect to 
save the event. 

A 250 keV threshold is imposed on each elec tron. T wo- 

are 



electron events from the 207 Bi deposit (Sec. IV Bl) 
eliminated by excluding electrons within a 4 wire radius 
of the source center. 

Occasionally, the fitter will have trouble converging to 
a good fit because of small defects in a visually acceptable 
track. The scanning software allows the user to remove 
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2 2.5 3 3.5 
Singles Energy (MeV) 

FIG. 14. Histograms of the 100 Mo sum- and singles-energy 
spectra for 2e~ events with tagged background decays re- 
moved. Events represented have single-electron energies 
greater than 250 keV and originate at least four wires from 
the 207 Bi calibration deposit. Also displayed are the results of 
the Extended Maximum Likelihood fit where the curves de- 
note total background (dotted), Pfov (dashed) and their sum 
(solid). 



these defects. This is done only when necessary, and 
always to the smallest possible extent. In the case of gaps 
in the track, a single added point will usually guide the 
fitter to the correct result. Ionization electrons formed 
where the track enters the 5 mm dead space between the 
TPC wall and the first plane of wires can arrive at the 
anode a few /is later, causing stray hits which may be 
safely deleted. An electron which scattered abruptly in 
the gas is usually handled by deleting an entire section of 
the track, leaving a single sinusoid for the fitter to work 
on. A status flag is recorded with the fit parameters 
indicating those few cases where the track was modified. 



C. Background Identification 



The tagging methods described in Sec. VI were used 



to eliminate, event-by-event, much of the U- and Th- 
series background, and estimate the number of back- 
ground events remaining among the /3/3 candidates. In 
the 100 Mo data set 16 background events were found, and 
in the 150 Nd data set 79 events, which otherwise would 
have been accepted as /3/3-decay candidates. The energy 
spectra of remaining events are shown in Figs. O and [l^. 
There are 410 remaining 2e~ events in the l0D Mo data 
set while the corresponding number for 150 Nd is 476. Ta- 
ble [y| shows details of how many background counts were 
identified and also summarizes the expected contamina- 
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FIG. 15. The 150 Nd counterpart of Fig. |l4| 

tion remaining in the data. The spectra still include the 
untagged background events listed in the last column of 
the table. 

The expected background contamination for 150 Nd can 
also be estimated from the activity derived from the lone 
electron fits. GEANT a-particle escape probabilities for 
212 Po and 214 Po (Table |) and 2e~ production coefficients 
(Table |T!l| ) were used in this calculation. The results, 
listed in Table [vj are in rough agreement with those in 
Table [v], with the exception of the Thorium daughters 
where the fit is less solid due to the interplay of several 
constituents in the same energy region. These fit results 
are not used in our final background estimates, but serve 
as a rough consistency check of the more precise deter- 
minations in Table 

In addition to U and Th backgrounds, isotopes outside 
these series contribute 2e~ events through the (3+ Moller 
process. Moller scattering by the entire lone electron 
spectrum can be estimated by Monte Carlo. Since a sub- 
set of these Moller events has already been accounted for 
in the contributions of isotopes listed in Table [y], those 
contributions are subtracted to yield the balance labeled 
"Moller" in the table. 

The method used to estimate the Moller production 
by the lone electron spectrum is as follows. The true 
number of lone electrons N(K), as a function of true 
energy K, can be obtained from the measured number 
N'(K') where K' is the measured electron energy. This 
is given by 



N(K) 



g(K) 



(i) 



where 



N' hiit (K)= / N'(K')P(K,K')dK>, 
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9(K) = 



N(K) 



j P(K, K 1 ) dK' = 1 is a probability distribution function 
while N^f ft (K) is the number of events surviving; both 
were determined from Monte Carlo using the measured 
100 Mo lone energy spectrum as input. The first iteration 
true spectrum, obtained in this manner, was input to the 
Monte Carlo to obtain the lone and 2e~ energy spectra. 
The former was found to be in excellent agreement with 
the measured spectrum; thus further iterations were un- 
necessary. 

The total (3 + Moller contribution, estimated from 
the Monte Carlo, is 21.8 ± 2.7 ± 2.3 for 100 Mo and 
13.9 ± 2.6 ± 1.5 for 150 Nd. The /3 + Moller events, re- 
maining in the data and listed in Table |v|, were obtained 
by subtracting the contribution of each isotope in the 
table. This contribution was determined by scaling the 
remaining events by the Monte Carlo determined fraction 
of j3 + Moller to total 2e~ events. 



D. Determination of /3/3-decay Rate 

The spectra of /3/3-decay candidates remaining are 
composed of real /3/3-decay events and a residue of 2e~ 
backgrounds. Good half-life determinations depend on 
accurate appraisals of the signals contained in the spec- 
tra. This can be done in several ways. 



1. Cuts Applied to Observed Spectra 

The majority of 2e~ backgrounds occur at low ener- 
gies, and can be greatly reduced by a high singles-energy 
threshold. A Monte Carlo simulation can then be used to 
estimate the number of events which would survive the 
cut, allowing us to estimate the number of events in the 
entire spectrum based on those events over our threshold. 

A 500 keV singles-energy cut is a reasonable thresh- 
old choice, as it eliminates a large portion of the back- 
grounds while retaining much of the /3/3-decay. Higher 
energy cuts would remove somewhat more background, 
but they rapidly degrade the statistics of the measured 
/3/3-decay spectra. This cut results in estimates for the 
standard uncut spectrum, of 420 ±34 events in the 100 Mo 
energy spectrum, and 480 ± 33 counts in the 150 Nd spec- 
trum. (This same cut applied to the | cos?? < 0.8 spec- 
trum implies 407 ± 42 100 Mo events, and 423 ± 39 150 Nd 
events.) 



2. Fits to the Observed Energy Spectra 

In addition to the above method, the data were fitted 
with the spectra expected from both the signal and the 
common backgrounds. This method allows full use of all 
data. 



The method employed for the fits was the Extended 
Maximum Likelihood (EML) method |l§j20) . EML anal- 
ysis uses one or more probability distributions which 
should describe the data in question, each of which has 
one or more free parameters describing its shape. The 
parameters are adjusted so that the probability of the 
data being derived from that distribution is maximized. 

For spectral fits, Monte-Carlo-generated (3(3 and back- 
ground spectra are used as probability distributions, with 
the only free parameters being their absolute normaliza- 
tions. Exceptions were 208 T1, 212 Pb and 212 Bi for which 
the measured 220 Rn injected spectrum was used. Open- 
ing angle inform ation was not utilized in the fitting pro- 
cedure (see Sec. [II B). 

The likelihood function is given by 



-AT 



(2) 



where P(xi,a) is the probability of observing event xt, 
given the distribution described by the parameters a, and 
Af is the normalization of the model. More specifically, 
the probability P(x, a) is the sum of all the probabilities 
of the normalized spectra of interest, i.e. 



where 



P(x,a) = ^2ajPj(x), 



Pj(x) dx = 1. 



The Pi(x) may include models of a variety of energy 
spectra, i.e. /3/32j,-decay sum, /3/32jy-decay singles, PPou,x~ 
decay sum, /3/3oi/, x -decay singles, 214 Pb sum, 214 Pb sin- 
gles, etc., as well as other distributions such as that of 
/3/3-decay opening angles. The parameters a are then ex- 
plicitly the absolute normalizations of the various model 
components, so that 



/v = £< 



When both singles- and sum-energy spectra are used in a 
fit, a single normalization is used to describe both spec- 
tra; i.e., there is a single parameter used to describe the 
number of events attributed to each process. With this 
scheme, the sum-energy spectrum normalization for a 
given background process would be a,j, while the nor- 
malization of the associated singles-spectrum would be 
2a j. 

A fitter was developed for maximizing Eq. (||) based 
on the CERN Program Library's Minuit minimization 
package. The fitter was designed for maximum flexibil- 
ity, and allows the user to select arbitrary combinations 
of models, apply various cuts to the data, restrict fits to 
specified energy ranges, and even apply a random smear- 
ing of the data for studies of systematic effects. 
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3. Inclusion of Background Estimates 

One has to choose a reasonable method for applying 
the various measured and calculated background models. 
Two extreme possibilities are either to fix all background 
model normalizations at the best estimates while leaving 
only the /3/3-decay model free to fluctuate in the fit, or to 
leave all models completely unconstrained. 

Neither of these extremes is particularly appealing. 
While the background estimates are based on careful 
measurements, and there is no particular reason to doubt 
them, the (3(3 candidate spectrum is, in a sense, an inde- 
pendent measurement of these contamination levels. Fix- 
ing the background models at some a priori normalization 
results in total disregard of the background information 
contained in the measured 2e~ spectrum. Similarly, al- 
lowing the model normalizations to vary freely ignores 
the best-estimates. In addition, the various background 
spectra can be qualitatively quite similar (e.g. the 214 Bi 
2e~ sum-energy spectrum and the theoretical /3/3-decay.) 
These similarities generate correlated parameters in the 
fitting procedure, adding non-physical constraints to the 
problem. 

A compromise to these extremes is to suggest the best- 
estimate values for the various backgrounds to the fitter 
as initial values, and then let the fitter vary the param- 
eters by some amount determined by the uncertainties 
associated with each. This is relatively easy to do in the 
Extended Maximum Likelihood scheme, since the energy 
spectra and the background estimations are independent 
experiments. This means that the corresponding proba- 
bilities may be simply appended to the likelihood func- 
tion. Specifically, the likelihood function is modified to 
include the probabilities of the background models tak- 
ing on given normalizations as determined by Gaussian 
probabilities centered on the best-estimates, with stan- 
dard deviations equal to the measured uncertainties: 



C = 



-(aj-a 3 ) 2 /2<r| 



where the aj are the best estimates for the normaliza- 
tions of the models, the tjj are the uncertainties in those 
estimates, and the Oj are, as before, the fitted parameters 
representing the normalizations of the models. 

The estimates summarized in Table ^ were used in 
an EML fit to the 2e~ spectra using the procedure de- 
scribed above. The uncertainties in the best estimates 
were scaled up from the l-cr values quoted in Table |v| to 
reflect 90% confidence levels. The only unrestricted pa- 
rameters were those associated with the ftfyv model nor- 
malization and the parameters of backgrounds not listed 
in Table The fits were performed on the standard 
spectra (i.e. no cosine or energy cuts), and consisted of 
a fit to all events in the sum-energy and singles-energy 
spectra. Table VII summarizes the inputs and results of 



TABLE VII. Summary of inputs and results of Extended 
Maximum Likelihood fits to the /3/3-decay 2e~ energy spectra. 
All uncertainties and limits represent 90% confidence levels. 



/3/3-decay Spectral Constrained Initial 



Fitted 



Isotope 


Model 


by £ 


Estimate 


Value a 


100 a it 

Mo 




no 






376.7 


+ 20.3 
— 22.4 




Moller 


yes 


20.5 ± 


5.9 


19.9 


+ 3.9 
— 3.9 




234 T-t 

A Pa 


no 






< 15.6 






Fb 


yes 


0.2 ± 


3.4 


< 2.3 






214r>- 

Bi 


yes 


0.12 ± 


0.17 


0.11 


+ 0.68 
— 0.11 




228 a 

Ac 


yes 


8.8 ± 


5.8 


8.5 


+ 3.9 
— 3.9 




208rjijb 


yes 


6.5 ± 


4.4 


5.9 


+ 3.4 
— 3.4 




211 D u 


no 






< 8.2 






Co 


no 






< 6.3 




150 at j 

Nd 




no 






414.4 


+ 22.3 
— 22.4 




Moller 


yes 


9.2 ± 


5.1 


8.9 


+ 3.7 
— 3.7 




d Pa 


no 






< 13.0 






214p b 


yes 


33.3 ± 16.7 


26.8 


+ 5.8 
-5.7 




214 Bi 


yes 


9.0 ± 


7.7 


9.4 


+ 4.5 
-4.5 




228 Ac 


yes 


11.9 ± 


7.2 


11.2 


+ 4.4 
-4.4 




208rpib 


yes 


8.7 ± 


5.5 


9.0 


+ 3.8 
-3.8 




211p b 


no 






< 5.1 






60 Co 


no 






< 3.6 






152 Eu 


no 






< 5.1 





the fitting procedure. 



a See discussion of uncertainties in Sec. |VIID3 
includes 212 Pb and 212 Bi. 



The uncertainties in the fitted normalizations are cal- 
culated by the MINOS routine within the Minuit min- 
imization package, and are derived by considering the 
shape of the likelihood function near the maximum. This 
procedure underesti mate s the uncertainties of the fitted 
components in Table VII because of the use of both sum 
and singles spectra in the fits. MINOS thinks it is work- 
ing with N + 2N = 3N independent events, but the 2N 
electrons from the singles spectrum provide no additional 
information for the overall normalization. We re-estimate 
the uncertainty in the number of (3(3 events by repeating 
the fit on the sum spectrum only. The resulting (3(3 com- 
ponent in this case is similar, so the sum-spectrum uncer- 
tainties for Npp are used in Table VIII and the half-life 
calculations. 

A second fit was performed on each data sample by 
removing models for which only an upper limit was ob- 
tained by the first fit. The fitted normalizations of the 
other models were found to be unaffected by this change. 

Figures |lj and [lB] show the fitted model spectra su- 



1G 



peri mposed on the experimental data. As described in 



Sec. VTID2, the fit considers each event independently; 



the binning shown in Figs. |lj and [l5] is for illustrative 
purposes only. 



4- Systematic Uncertainties in EML Fit 

The EML fit is a point-by-point fit, rather than a fit 
to a binned histogram; this feature makes it simple to 
estimate the effect of the TPC's energy resolution on the 
fitted results. The energy of every electron in the data 
set has an uncertainty associated with it, so by randomly 
altering each of these energies in keeping with a Gaus- 
sian distribution defined by the best-estimate energy and 
its uncertainty, a hypothetical energy spectrum can be 
generated from the altered experimental data. This new 
spectrum can then be analyzed with the EML fitter, and 
will probably produce results slightly different from the 
unaltered data; these differences are due solely to the en- 
ergy resolution of the experiment, not the statistics of 
the data set. This procedure was repeated about 1000 
times on both the 100 Mo and 150 Nd data sets, and dis- 
tributions of /3/?2i/ fits were generated. The widths of the 
distributions were quite small, and imply an uncertainty 
of only about 0.5%. 



TABLE VIII. Parameters used in half-life determination. 
The "Source mass" is the Mo or Nd2 03 material actually 
used, not just the /3/3-decay atoms. "Observed /3/3-decay 
events" is based on maximum likelihood fit to data; see text 
for description of other methods and results. Uncertainties 



are la. 


1 cLI dill 




Parameter Value 


ctcr 


Description 


100 Mo Run 


150 Nd Run 


r~ a 


2e~ efficiency (%) 


10.1 ±0.1 


11.5 ± 0.1 


fent 


Isotope enrichment 


0.974 


0.91 


/mol 


(3(3 atoms/molecule 


1 


2 


Mtot 


Source mass (g) 


16.7 ±0.1 


15.5 ±0.1 


w moX 


Mol. wt. (g/mol) 


99.91 


346.4 


N A 


Avogadro # (mol -1 ) 


6.022 x 10 23 




Live-time (h) 


3275 ± 2 


6287 ± 3 


N?f3 


Observed (3(3 events 


377±2g ± 2 


414±25 ± 2 



TO % systematic uncertainty not included. 



E. Half-Life Calculations 

The half-life of the /3/3-decay isotope under investiga- 
tion is given by the fundamental quantities determined 
in this experiment as 



12 ~ (loo) /( 



f ( Mtot \n ( tlive 

enr Jmo\ TTZ JVA — 

\Wnxo\J \Nf3fj 



In 2. (3) 



The parameters in Eq. (|J) are described in Table VIII , 
which also summarizes their values. Using these values, 
the EML-based half-lives for the two isotopes studied are 



rr2iy /100 



l/2\ 
r i/2( 



Mo) 



(6.82+^* ± 0.68) x 10 18 y, 

s 

y- 



Ti% ( 150 Nd) 



(6.75lo.42 ± 0.68) x 10 18 



A comparison with other experiments can be found in 
Table and with theory, in Table [x|. 

The estimated number of /3/3-decay events based on 
single-electron energy threshold cuts may also be used in 
lieu of the EML-based estimates in the half-life calcula- 
tions. The half- lives derived from those values are plotted 
in Fig. 16. The hatched regions in the figure represent 
the result of the EML half-life calculation. The EML 
result agrees nicely with the the singles-energy thresh- 
old method near the favored threshold of about 500 keV. 
This figure does not include systematic errors, which con- 
tribute approximately the same uncertainty to each mea- 
surement. 



TABLE IX. Comparison of the results produced by this 
work to those from other experiments. Uncertainties and 
half-life limits are at the 90% (68%) confidence levels. Limits 
and uncertainties in square braces indicate that the confidence 
level was not specified. 



1 1/2 

(10 18 y) 


1 l/2 

(10 21 y) 


rp0u,x 
1 l/2 

(10 20 y) 


Group 


100 Mo Run 


6.82+ ± (0.68) 
11 6 +(3 - 4) 

Xi - D -(0.8) 


> 1.23 

> 2.16 


> 3.31 

> 3.9 


This work 

uci p 

UCI 


r +5.4(3.0) 
XJ -° -2.8(2.0) 
o o + (2.0) 
°° -(1.0) 

10 ±[0.8] ± [2] 


> 2.6 

> (0.71) 




Osaka HQ 
INR H 
NEMO H 


9.5 ±[0.4] ± [0.9] 


> 6.4 


> 5 


NEMO 2 H 


150 Nd Run 


6.751 £$±(0.68) 


> 1.22 

> 2.1 


> 2.82 

> 5.3 


This work 
UCI 


17 1 5° ± 3.5 






ITEP/INR @ 
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TABLE X. Theoretical half-lives for 100 Mo and 150 Nd. 
Zero- neutrino predictions assume < ra„ > = 1 eV. 



Isotope 


1 1/2 (y) 


1 1/2 (y) 


Ket. 


100 Mr. 


i i ^ v i n 18 

I.IO A 1U 


1 97 v 1 n 24 

l.Zi A 1U 






2.87-7.66 x 10 18 


2.6-4.7 x 10 23 


@ 




6.0 x 10 18 


1.9 x 10 24 


§] 


150 Nd 


7.37 x 10 18 


3.37 x 10 22 


§] 






1.05 x 10 24 


u 




6.73 x 10 18 




§] 



F. Kurie Plots 

The Kurie plots were generated by first subtracting 
the background 2e~ sum energy spectra, from the candi- 
date spectrum, using the amounts given by the EML fit 
results. The next step involved obtaining the true 2e~ 
residual spectrum. Here Eq. (|l|) was used where K and 
K' now refer to the true and measured sum energy and 

9{K) = IukV' 

where N r (K) is the true number of events in the 
sum energy spectrum with singles energies greater than 
0.25 MeV. Also, for this calculation, P(K,K') and 
N^ it (K) were determined using the theoretical (3(5 spec- 
trum as input to the Monte Carlo. (This was later re- 
placed by a randomly generated spectrum in order to 
determine the systematic uncertainty.) Finally the Kurie 
plot values were calculated from the true residual spec- 
trum using the formula [ p8| 

(Q/3/3 - K) oc ((dN/dK)/{(K - 2r)[h(K) + f T (K)}})^ 5 
where 

fo{K) = K 4 /30 + K 3 m/3 + 4K 2 m 2 /3 + 2Km 3 + to 4 , 
f T (K) = t(K - t)[K 2 /15 + 2Km/Z + 2m 2 /3 
+ t(K-t)/5], 

and to is the electron mass. This formula was ob- 
tained by including the singles threshold, r, and using 
the Primakoff- Rosen approximation for the Coulomb ef- 
fect; the approximation produces < 1% distortion over 
the plotted range. The Kurie plot (see Fig. [l?]) projects 
to an end-point energy for the 100 Mo source of (3.03 ± 
0.02 ± 0.06) MeV, in good agreement with the mass dif- 
ference of 3.03 MeV. The corresponding 150 Nd value is 
(3.44±0.02±0.06) MeV, also in agreement with the mass 
difference of 3.37 MeV. 
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FIG. 16. Graphical representation of half-lives derived 
from the number of events at various singles energy thresh- 
olds. (Open circles are shifted by -0.01 MeV for clarity.) The 
hatched band represents the result of the the EML fit to the 
energy spectra. Uncertainties are statistical, and represent 
90% confidence levels; systematic uncertainties contribute ap- 
proximately 10% in addition to those shown. (One should 
keep in mind that nearby points are strongly correlated.) 




Energy (K) (MeV) 




0.5 1 1.5 2 2.5 3 3.5 



Energy (K) (MeV) 

FIG. 17. The Kurie plots, for 100 Mo and 150 Nd. Also 
shown are linear fits to the points and the x-intercept indicat- 
ing the end-point energy. (F(K,t) = (K -2T)[f (K) + f T (K)]. 
See text.) 
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G. Exotic Decay Limits 

1. PPov -decay Half-Life Limits 

The expected Pftov distributions were generated with 
a Monte Carlo simulation. For 100 Mo, the FWHM lies 
between 2.89 and 3.18 MeV and for 150 Nd, this region is 
bounded by 3.21 and 3.54 MeV. There were no counts 
from either isotope observed within the cut regions, and 
poisson statistics allows us to reject 2.3 events at the 90% 
confidence level. 

Half-life limits are derived from the above limits by 
ap plying Eq. (ph. Two of the parameters listed in Ta- 
ble VIII must be altered for this calculation: The 90% CL 



decay limits are used in place of the number of observed 
decays, and a new efficiency for the 0(3q v case is given 
by Monte Carlo as epp = 11.12 ± 0.07 for 100 Mo and 
epp = 11.50 ± 0.07 for 150 Nd. With these parameters, 
the half-life limits are (at 90% CL) 

T%( 100 Mo) > 1.23 x 10 21 y, 
T 1 ° / I ' 2 ( 150 Nd) > 1.22 x 10 21 y. 

(See Table |x| for theoretical predictions.) These /30o„ 
half-lives do not give competitive limits on neutrino mass 
and right-handed current parameters. 



this case, negative normalizations are non-physical solu- 
tions) , so a direct integration of the likelihood function 
must be performed in order to calculate the confidence 
intervals. 

The 90% confidence limit on the number of /3/3o^,x 
events in our data sets is 7.76 events in the 100 Mo data 
set and 9.73 events in the 150 Nd data set. Half-life limits 
can be calculated by using Eq. (^) with these decay limits 
in place of an observed number of decays, and using the 
appropriate /3/?oi/,x efficiencies of epp — 10.13 ± 0.06 for 
100 Mo and epp = 11.25 ± 0.07 for 150 Nd, 



1, giving 



T y,x ( ioo Mo) > 3 31 x 1Q 2o Yj 



T 



^ x { 150 Nd) > 2.82 x 10 20 y, 



at the 90% confidence limit. The (3(3q v , x half-life limit 
published by our group in 1994 Q was based on an anal- 
ysis of binned data, and the fitting routine has since im- 
proved. We consider the current analysis to be more reli- 
able. We use the matrix elements of Ref. j29] to put limits 
on the effective neutrino-majoron coupling constant and 
obtain 

<g v , x > ( 100 Mo) < 6.26 x 10~ 4 , 
<g„, x > ( 150 Nd) < 9.96 x 10" 5 , 

at the 90% confidence level. 



2. /3/3ov, x -decay Half-Life Limits 



The expected energy spectra for the 0ftou,x mode 



are continuous distributions, spanning the entire energy 
range available to the Pfcv channel. As such, lim- 
its on the contribution from /3/?oi/,x cannot be deter- 
mined by simply counting events in a predetermined re- 
gion, and the maximum-likelihood method described in 
Sec. VII D 2 must be used. We restrict our analysis to 
the "ordinary majoron", index=l in ref. ||. 

Maximum likelihood fits were performed on both the 
100 Mo and 150 Nd data sets, using the same fitting pro- 
gram previously used for the j3p2v analysis. In this case, 
the fit was restricted to energies well above the majority 
of /3/?2i/ backgrounds, where the majoron decay signature 
is significant (i.e. singles-energies greater than 750 keV, 
and sum-energies greater than 1.5 MeV). The fit was 
performed jointly over both the sum- and singles-energy 
spectra, and included only two models: the fifiiv and 
Pf3ov,x Monte Carlo energy spectra. Under the assump- 
tion that poor energy resolution at high energies will ad- 
versely effect the reliability of our /3/3a VjX limit estimates, 
we cut the data at | cos#| <0.8. 

In both the 100 Mo and 150 Nd cosine-cut spectra, the 
EML fitting procedure indicates that the data sample 
is consistent with there being no (3Poi>,x events in either 
spectrum. The algorithm used by Minuit to determine 
confidence intervals, MINOS, cannot be used when the 
most likely solution is against a physical boundary (in 



VIII. CONCLUSIONS 

Although our choice of a very thin source plane pre- 
cluded sufficient source mass to give competitive limits 
on P/3q v , it did result in exceptionally clean Pfcv energy 
spectra composed overwhelmingly of /3/3-decay events. 
These spectra were produced without subtraction of any 
untagged background events. The derived Kurie plots 
are straight over a broad energy range, and intercept the 
energy axis close to the expected values. 

Improvements in fitting techniques and in determina- 
tion of the detector efficiency have resulted in a 100 Mo 
half-life that is shorter than we reported previously [Bll . 
We now believe that the maximum-likelihood fit per- 
formed on that early data set significantly over-estimated 
the background contamination, contributing to an artifi- 
cially high half-life estimate. 

The intriguing high-energy anomalies in our earlier 
/3/3-decay sum-energy spectra did not withstand the im- 
proved energy resolution that resulted from a doubled 
magnetic field strength. No suggestion of f3(3ov. x remains 
in either sum spectrum. At our level of sensitivity, the 
double beta decay phenomenon is well described by stan- 
dard theory. 
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